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To remind in brief a discussion which is already distributed over several publications: 

A damping mechanism of giant resonances (GR) is well established and represents now basic knowledge in nuclear 
structure physics. Calculations performed by many groups of authors within different microscopic approaches confirm 
that a spreading width (due to a coupling of collective modes, phonons, to complex configurations) is the main part 
of the total GR width in medium and heavy nuclei. In light nuclei, a coupling to continuum (an escape width) also 
■ plays an essential role. 

A damping mechanism of GRs in a phenomenological phonon damping model (PDM) in its PDM-1 version is 
different from that (see, an important clarification in ||). A collective phonon fragments within PDM-1 as a result of 
ON ' coupling to simple and not to complex configurations, i.e. only the so-called Landau damping mechanism is accounted 
for. A coupling strength is a phenomenological model parameter which is adjusted to reproduce the GR width known 
from experiment. Agreement with data provided by fits within the PDM may be defined as from very good to 
excellent. 

In a recent article [[ll| which raised the present discussion, it has been concluded that these type of fits confirm "the 
ly*-) [ authentic damping mechanism" of the PDM as "the result of coupling between collective phonon and non-collective 
p-h configurations" (i.e. the well established knowledge on the GR properties was put in doubt). This conclusion 
\ has been criticized in my article |Q. It has been argued that this model has the Breit-Wigner (BW) form for the 
phonon distribution as an ad hoc input and thus, even excellent description of the data available is not surprising. 
A fruitfulness of an idea to make conclusions from fits in which model parameters are adjusted to described physical 
, observables has been put in doubts. 
& ' Although my evaluation of the PDM in |1 was made for the point of view of general physical grounds, Dang et al. 
did not agree with me in the forthcoming publication [0] . They claim that I consider some specific case ( "his case" ) 
which cannot be attached to the PDM and all my arguments "are either wrong or irrelevant' . I cannot agree with 
^ ' their conclusion and present below additional arguments in a sequence following the paragraphs in || : 
£ ■ 

2. For the giant dipole resonance (GDR), the energy scale associated with variations in a coupling matrix between 
. £h ' a phonon and uncorrelated Iplh states is of the order of a few hundred keV. The width of the GDR strength function 
is of the order of a few MeV. So, I do not agree that the condition cited in [|| from || is satisfied in the GDR region: 
' why are a few MeV small compared to a few hundred keV? 

I know only one PDM-1 article jl2| in which it is assumed that a phonon interacts 40 times stronger with some 
specific configurations than with other ones (see more on this article in 9. below). In all other PDM-1 papers, we 
find a single phonon which interacts equally with all Iplh configurations. I do not want to discuss here the PDM fits 
at non-zero temperature. To keep on reproducing the data in hot nuclear, Dang et al. have to assume for unclear 
reasons that a phonon prefers to interact with lplp and \hlh configurations about 10 times stronger than with Iplh 
configurations. Again, as in the case of cold nuclei, an idea to provide the best fits is preferred to understanding of 
the physics. I think it is a blind way for theory. 

It is true that PDM equations are presented in a general form in many papers by this group with different V qiSl . 
But the point is that they are never used in actual calculations in this form. For this reason, I prefer to discuss what 
is used in calculations rather than what is written and not used even by the PDM authors themselves. 



* I cannot not neither confirm nor disprove that an article |l| was rejected by PRC as a reader learns from [pj. It is common practice that 
a manuscript submission to any scientific journal remains a confidential information of the author(s) until the manuscript is accepted 
for publication. In PRC, it is protected by a secret accession code. It seems to me at least not ethical to publish openly on Web a 
confidential information of another person(s) as Mr. N. Dinh Dang did in foj . An editorial decision may be published only with a formal 
permission from the Editors. I doubt that Mr. N. Dinh Dang has it. I also interpret this fact as an attempt to influence a reader by 
non-scientific arguments in a scientific discussion. 
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3. It is very simple to transform Eq. (1) in pi for m)^ into Eq. (1) in |jj for W2, although Dang et al. claim it is 
impossible. For that, one needs to switch off an additional PDM smearing, i.e., consider the limit e — > 0. This would 
bring immediately to the first line of Eq. (2D-14) in ||. Eq. (I) in (for a constant coupling strength) or its general 
form in ||: 

w 2 = Y^ty™) 2 ( 2D - 14 ) 

a, a 

for the second moment W2 is relevant to the PDM as well as to any model which deals with interactive systems. 
Of course, to perform this transformation one should use the PDM strength function introduced in Ref. fl3]| : 

S q {E) = I (I) 

*{E-u> q -P q {E)) 2 +1 2 q {E) 

where r y q (E) is the PDM damping, P q (E) is the polarization operator (see, e.g., Ref. for definitions), and uj q 
is a phonon energy, a model parameter. The strength function S q {E) presents fragmentation properties of a PDM 
phonon over eigen-states of the PDM Hamiltonian smeared with an additional parameter e. Parameter e appears in 
5(E) = e/[n ■ (E 2 + e 2 )} for ^-functions in ~f q (E). 

I point this out because the st reng th function (|l|) has been replaced in the forthcoming PDM articles [D], [l2| [Tj|, 
|l|, ||, |l|, m, |(], || ||, |4], by its approximate form: 

S>(E)= l - ^ (2) 

v{E-E GDR f +1 2 q (E) 

where Eqdr should be taken as a solution of 

f{E) =E-tu q -P q (E) = . (3) 

Eq. (||) has been obtained from Eq. (0) by expanding P q (E) near a solution of Eq. (^), Egdr, and then extrapolating 
the properties of this approximation far away from Eqdr.- In the limit e — > 0, Eq. (||) has 7V+1 solutions corresponding 
to eigen-energies of the PDM Hamiltonian. 

4. I never claimed that the BW form for the phonon distribution is assumed within the PDM. But it is indeed an ad 
hoc input for PDM calculations. I may refer again to || where we read that "the Breit-Wigner form for the strength 
function is an immediate consequence of the assumption of a constant coupling to the other degrees of freedom of the 
system" . The BW under discussion has nothing to do with definition of the PDM strength function. Indeed, in the 
limit e — > 0, S q (E) turns into a set of infinitely narrow lines while their envelope still remains the BW. 

5. I do not agree that the calculation with random values of E a in |IJ "no longer corresponds to the PDM". I 
have used the PDM Hamiltonian and details on a spectrum and model parameters are only technical details of a 
calculation. The purpose of my calculation is to demonstrate that 11 the crucial feature of the PDM is the use of 
realistic single-particle energies" is of marginal importance when a configuration space is not small; everything is 
determined by the BW discussed above. 

Eq in Ref. ^belongs to the Lorentz line in a hypothetical nucleus and not to my PDM fits. Eigen energies in my 
calculation in Ref. [Q were obtained from Eq. (||) in the limit e — > 0. 

6. I agree that if something "is by no mean[s] obvious" it has to be checked. My experience of microscopic 
calculations tells me that the increase of collectivity tends to the increase of a coupling strength. Of course, it is not 
necessary that everybody should trust my experience. But then, there are no other alternatives: the one, who puts it 
in doubt, should check it independently. 

7. I never claimed that there are some reasons "why the values of fx for 40 Ca and 48 Ca should be the same" as 
there are no reasons to keep this parameter fixed along chains of isotopes. As pointed out in Q, this parameter has 
no physical meaning. My issue is that one cannot learn anything from agreement with experiment from fits in which 
a free parameter is adjusted to a described observable. It is important to stress once again: 



The Phonon Damping Model (PDM) has three phenomenological parameters, u> q , fx, and ci for 
position, width, and amplitude of the GDR, respectively, and the Breit-Wigner as an ad hoc input 
for its shape. Numerical values of these parameters cannot be determined from independent data. 



On the page 4 of the article under discussion jfl]], we find: "For double closed-shell nuclei 16 O and 40,48 Ca, where 
the pairing gap is zero, such a kind of enlargement of configuration space is compensated simply by a renormalization 
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o//i, which reduces its value by ~ 25% for 16 0, and ~ 35 — 37% for 40,48 Ca." I read this statement after discussion 
of open-shell nuclei as a renormalization of fi in closed-shell nuclei in respect to open-shell nuclei. Or do Dang et at 
mean to say that calculations in double-magic nuclei have been performed with pairing as one may conclude now from 
P|: "The results for GDR in 16 have been obtained already within the enlarged space"! Obviously, only the authors 
know whether they renormalizcd their f\ in calculations along chains or not. Of course, I will take out the statement 
on the /i renormalization from 16 to 18 if it is not true. But before that, I need some help from the authors as to 
how a reader should interpret the above cited statements. 

8. It is clearly explained in Q why comparing the PDM predictions in 40 Ca to the data ||, the strongest 1~ state 
observed should be excluded from consideration (because it has a two-phonon nature and two-phonon states are not 
included in the PDM model space). Thus, the PDM 0.25% of the TRK EWS corresponds to 0.007% and not to 
0.025% from this experiment. It seems to me that Dang et al. try to hide again a huge disagreement by misleading 
comparison. 

The same conclusion, that the PDM is not capable to reproduce "the significant experimental difference in the El 
strengths" 40 Ca/ 48 Ca, has been obtained independently by another group of authors 0. As they write, "It is important 
to note that the parameters of the PDM are adjusted to reproduce the gross structure of the GDR while investigations 
of ^-ray strength function models show that the extrapolation of the strength distribution down to energies below the 
particle threshold leads to unrealistic high dipole strengths and overestimates the experimental data" . Thus, it is not 
only my point of view that a conclusion in |ll| (on a quantitative description of pygmy dipole resonance within the 
PDM) is not justified. 

It is not true that the PDM with a structureless phonon has no problems with double counting. If a phonon internal 
structure is not accounted for, it does not change the physical meaning of the phonon. The PDM configuration space 
contains a phonon and uncorrelated lplh configurations. The last ones are also excited from the ground state and each 
of them has its own B\ p ih(El) value. If lplh spectrum is rather complete (it is always true in the PDM calculations) 
these uncorrelated lplh state alone exhaust about 100% of the TRK EWSR. But the PDM physics is determined only 
by the phonon strength function and not by its sum with N strength functions of uncorrelated lplh configurations. 
This is equivalent to Bi p ih{El) = within the PDM. 

9. The previous article on pygmy resonances by Dang et al. ]l2] | was not a subject of (TJ. But if Dang et al. raise 
discussion on it in 0,1 have some comments on it too: 

The capability of the phonon damping model (PDM) to describe giant dipole resonance (GDR) damping in neutron- 
rich nuclei has been tested in [fl2"[ . To mimic essential differences of double-magic and exotic nuclei, a coupling between 
a phonon and some lplh configurations near the Fermi surface has been strongly enhanced in the last ones. As a 
result, a phonon interacts with these selected lplh configurations with a strength "equal to 41 MeV for oxygen 
isotopes, 13.856 MeV for calcium isotopes, and 6.928 MeV for tin isotopes" M. Let us try to understand how it is 
possible to stand such an enormous coupling strength which is far away from nuclear structure scales and report an 
agreement with experiment from this type of calculations. For that, I have repeated the PDM calculations for 16 and 
18 at zero temperature keeping all the details of Ref. jl^] and employing realistic lplh spectrum from Hartree-Fock 
calculation with SGII Skyrme force s B . 

The results of my calculations fllCfl are shown in Fig. [|. A difference in S q (E) for 16 O and 18 O in Fig. |] is 
dramatic but not surprising. It is due to the fact that a phonon couples to all lplh configurations with an equal 
strength of F% = 1.025 MeV (a PDM parameter) in 16 while a coupling strength for [1^1/22^5/2] v configuration at 
E[i Pl/2 2d 5/2 ] l , — 8.2 MeV in ls O has been enhanced to F[ = 40i*\ = 41 MeV following the details of calculations in 
Ref. Q. As a result, we find the GDR in 18 O between 40 and 60 MeV and about 40% of its strength is pushed to 
—20 MeV. The energy of the ground state is MeV. 

Comparing dashed lines in Fig. |], one may be surprised that S' q {E) does not feel this enormous matrix element of 
41 MeV in 18 confirming the results in [JL2J. But a very important detail is that S' q (E) in Fig. [j] has been calculated 
with Eqdr = 22.5 MeV from Ref. |Q. It has been done in an attempt to reproduce the GDR strength functions 
published in Fig. 3 of |L2| which are found in agreement with the data available in 16 and 18 0. Taking into account 
that the employed lplh spectrum might be not exactly the same as in Ref. |l^], it is possible to conclude that dashed 
curves in Fig. [j] reproduce the results in Fig. 3a, 3b of Ref. |Q on a rather good qualitative level. 

The only problem is that it is not possible to obtain Eqdr = 22.5 MeV in 18 O, reported in Ref. |Q, with parameters 
from this article as a solution of Eq. (^). To demonstrate this, let us consider a behavior of the function f(E) of 
Eq. (|3|) in 16 and 18 0. In 16 0, it has a tendency of a continues increase with fluctuations reflecting lplh poles 
smeared by the parameter e and crosses y = line in my calculation at Eqdr = 18.5 MeV (see, Fig. ||, left). In 18 
(Fig. ||, right), a fluctuation around -E , [ip 1/2 2d 5/2 ]„ increases enormously because of 41 MeV coupling matrix element 
corresponding to this pole yielding a spurious solution of Eq. (||) at this energy. The physical PDM solutions of 
Eq. (|) in 18 have E GDR = -19.7, 50.6, and 55.0 MeV which are very different from E G dr = 22.5 MeV, reported 
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FIG. 1: Strength functions S q (E) (solid line) and S' q {E) (dashed line) of PDM phonon distribution in 16 and 18 0. 
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FIG. 2: Functions f(E) in 16 and 18 0. Vertical line is plotted at 22.5 MeV. 



in Ref . |T^] . S' q (E) calculated with any them is dramatically different from the one in Fig. | (right). 

The " 18 effect" can be obtained even without any calculations. For that, one may neglect a phonon coupling to 
all Iplh configurations (with a "weak" matrix element F±) except for [1^1/22^5/2]^ configuration. Then, in the limit 
e — > 0, Eq. (pi) transforms into quadratic equation: 



E 



E- Ei 



= 



[lpi /2 2d 



(4) 



5/2) 



where a factor n accounts partial occupation of v2d 5 / 2 level in 18 0. Eq. (Q) yields the PDM eigen states at —19.0 and 
49.4 MeV with a phonon strength distribution among them as 40% and 60%, respectively. 

It becomes clear that agreement with experiment for 18 O (and accordingly for other neutron-rich nuclei) reported 
in Ref. [Q has been obtained by making use of the approximate PDM strength function S' q (E) and the GDR energy 
which is not a solution of Eq. (g) as announced. Correct PDM strength function with parameters from Ref. |Q for 

is presented by solid curve in rig ht part of Fig. 0. 

1 and 10. I have examined the PDM from the point of view of general physical grounds. My arguments and 
conclusions are presented in jlj and above. I think a reader may independently conclude whether general rules are not 
for this model (as the claims of Danger al. in [H may be understood) and whether one learns any physics from the 
PDM fits in Refs. @ H 0' EH HI S H 13 H H H II 11 although agreement with experiment 
is always reported by the authors. 
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